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I. INTRODUCTION 

In pursuit of an interest in the problems and potential 
of sateilite tracking of earth's surface or low altitude 
targets, the interested student quickly discovecs a plethora 
of information available on the subject, but little in the 
way Of analysis of the entire problem, start-to-finish. fhe 
purpose of this report then is to provide thateandiycie, 
beginning with the mechanics of the orbiting satellite; next 
developing a model to represent the target notion; and 
continuing tnrough an observability analysis of tke prcklen. 
Finally, a first-cut attempt is made, based on the observ- 
ability study, to design a deterministic nonlinear observer 
and an Extended Kalman Filter (hereafter referred to as an 
EKF) to observe and track the target. 

Chapter II of this ceport examines in depth the physical 
dynamics of a satellite orbiting the Earth, and then the 
target moving on or just above its surface. A working nodel 
is developed for the satellite, incltding a foci 2c. 
ground track; its rotation about itS own axeS aS it roves 
through space; and the translation of its raw observed data 
into eafth's surface eoctds fave Roth position and rate 
data are considered, and specific formulae are deveiozed 
alonc the lines of a computer alyorithm for coordinate 
trans lormart on. 

Chapter Iii deals more specifically with the modei for 
tne target itself. Actually, two models are considered: 
first, a realistic model, necessarily nonlinear, due tc the 
hnon-synchronous moticn of the target and satelite me: 
concentric spheres; and then secondly a slightly nonreal- 
istic iinearized version. Both surface and alr targetouarc 
considered, Lut for brevity's sake, only the surface target 


is discussed fully in the variots examples. 
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Chapter IV finaliy tackles the problem of state observ- 
ability. Observability, as defined here and in [Ref. 11, is 
the capability, given the observations and knowledge of the 
control input to the system over tine, to reconstruct the 
entire system state kFack to the initial condition. Testing 
BOEw tis CONGHti1Oh On avlinear Ssystcn is perfectly 
BOE4araG; the method is fuily ceveioped in [Re=. 2] and eise- 
where. The opservability of the nonlinear systen 
the mere interesting problem, as the problem is mere reai- 
istic, andthe methodology for testing of nonlinear opserv- 
ability is far less well develcped and/or documented. ee) 
that end, the method which is aewly developed in [Ref. 1] is 
herein employed, and the analysis uitimately becomes as nuch 
of a test of the methcd as the cbservability of the systen. 

Based on the results of the analysis in Chapter IV, 
Chapter V is then an admittedly imperfect, first-cut attempt 
at designing a deterministic ncniinear observer, and then 
tracking the target with an EKF. The comments in [Ref. 3] 
concerning filter initiaiization and measurement covariance 
eStimation recome pertinent. 

Finally, Cnaapter Vi su@marizes tne thesis' rcesuits, and 
makes some suggestions for further aualysis, snecificaliy ia 
426 aréa oO: multzple taryetsS and onpservers. Pee TOA wo sre 


Steere cance OF the revort 1s also ciscussed. 


II. SATELLITE DYNAMICS AND OBSERVATION OF TARGETS 

In considering the tracking of a SUEFiaGe OERat pee oese: 
by a satellite kurtling through space at speeds approachiny 
15000 knots, it is obviously imperative to understand the 
dynamics of the moticn of the satellite. Only then can a 
relative coordinate system be established ;withmene sue ama 
goal of describing the target's motion over the Earth's 
surface; first, however , the satellite's own motion must be 
Similarly described. 

In the following derivations, the Farth will be approxi- 
mated py a uniformly-dense, perfect sphere. This approxima- 
tion has often been characterized cualitatively as being 
twice as good as that of a "sphericai’ basketball; its cguan- 
titative implications wiil be discussed in Sections B and C 


below. 


A. SATELLITE REFEREACE ORBIT 


Lagrangian Mechanics, as developed in [Ref. 4] offers 
the easiest method for anaiyzing the satellite's moticn. 
The Earth and the satellite can each be represented by pcint 
masses; in the absence of external disturbances, it can be 
shown that the sateilite's orbit will remain planar. the 
gravitationai potential ofr the Sateilite is inversely 
prorortional to its distance £Lecm the Satin. Usirgeer emt 
polar coordinates as shown in Figure 2.1, this potential 


energy V is then described by 
_ — GM / 
\V = C 


wnere G 1S Newton's constant of universal GEayrcatton eo 


) 


the mass of the Earth, and m 8S the tasSaseort themeetc ha. 
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Figure 2.1 Satellite Orbit in its Plane 


The circular orbit is the simplest case, Sicierenadt the 


velocity is purely tangential and can be characterized as 


\l, . coq 


the kinetic ee is then 


¥ “le 
1 ( <j & : *h 
| = aie << s ‘a gael ) 
and in the absence of any external disturbances, the 


iaAgmangian cunction is written as follows: 


3 a ee, 
SE THV= deere e SM Ve, 
By the rules of Lagrangian mechanics, 
d /oy\. of 
alex = Sie 7. 
or po lea 
aap 


pes, 


OO Gag chy 
or 
CO Ts 
AEM O] = 0. (on 
ae Wsleshoy aa Ye indicates’ that angular womentun is 


conserved (i.e. its derivative is zero), and that @ will be 
constant for a given orbital radius. FULth CE, =e 


angular momentum vector iS interpreted as 
—, - = zi 

) : on (Fx) Sn UF 

, 


then the fact of its conservation indicates that the orcit 
will remain stable in its original plane. The“ antEre oduct tog 
of a third coordinate, which ultimately becomes the angie of 
BB Gee Jest joYs ie aleyel © ak. and which remains arbitrary throughout the 
Lagrangian derivation, proves this conclusively. 

Solving equation 2.1 for the the physical constants 


reveals that 


Jao 
vo = GM. 
Considering that @ ia 


wnere T is the satellite's pericd, then 


3 
eo Gt 
Vo at * 
for a satellite in an earth orkit; this is Kepler's Seccnd 


Law Of OEbatal Motrvons 
The radius in this equation is the distance of the 
satellite from the Earth's center; for altitude arove the 


Earth, the Earth's radius, here represented as 


ie 3438nm = 6366kn 
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must be subtracted, i.e. 


Of special interest is the satellite at Zero altitude, 
aiso knewn aS a derget satellite, with a period of 84.5 
Minutes. As [Ref. 5] indicates, this gives a value for the 
lower bound of a satellite's period. Some other periods and 


Geearudes OL interest are summarized in Table I=: 


—_ °°» TS Di eS A 


TABLE I 


1 


Typical Periods, Aititudes, and Orbital Velocities 


| 
| Period T Altitude h Orbital Velocity 
| 
| 


SOP ain 150 ar 15930 knots 
hr ng nots 

2 23 OY) 11930 knot 
Omir 5600 no 9470 knots 
UZ DE 190900 nia HOU MRenO LS 
3 Ms 22 I/O Knots 


(Ces fii aie Gee Soo geen Ge Gee Gi eS Ge Bee 


n 
23enc 56min no 


iies@last orbalt, @with a period of one Solar day, is the 


| 
hn 


geosynchronous orbit frequently adopted in satellite conmmu- 


mreations apdlications; for tracking of surface-bound 
targets, however, low-altitude orbits will be of much 
greater interest. Typical altitudes will generaily be iess 


than 1000 nautical miles, depending on the sensitivity and 
accuracy of the tracking sensor: this corresponds to 
periods of two hours or less. 

Finally, it can Ete shown by the method of the potential 
nergy well [Ref. 6], and alsc by the application of the 
theory of perturbations [Ref. 4j, that this reference crbit 
will be stable in the presence of small disturbances in 


€ither the orbital plane or the third dimension, and these 
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disturbances will produce oscillations about the reference 


orbit with its same period, i.e. 





Be. SATELLITE GROUND TRACK 


The earth-sateliite model as developed in Section A is 
obviouslv over-simplified for the purposes of target 
tracking; thus far, the Earth has been represented only as a 
point mass rather than as a spinning globe, which would make 
it extremely difficult to describe any motion across’ the 
Earth's surface. 

One intrinsic feature of a satellite's orbit used in 
measuring its flight path is the inclination angle i, which 
represents the fositive angle between the orbital plane and 
the equatorial plane. This angle can be used, together wit! 
earth's surface coordinates, to describe the path over 
ground, or ground track, of the satellite. 

In Figure 2.2, the Z axis coincides with the Earth's 
polar axis; L and A measure necrth latitude and east lonji- 
tude respectively; andr, © , and vV, describe the satel- 
lite's motion in its orbital plane as before. 

Tf the satellite is considered to have originally 
orbited in the equatorial plane, its positior would have 


been described by 


A = © ces, 
Yr rong, 
=O 
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regure 2. 2 


mee the orbital 


through i degrees, 


thet is, now 





plane 


— 


ad 


Inclination of the Orbital Plane 


is tnen rotated about the X 


then the coordinates are transformed by 


| © & 
© CCS Fe as - 
) 
= Sig 4 OSC 
e ees) 


Sa d Ce, -; 


SIN a) Si (oe) 


and the latitude and longitude are described by 


i 





(es Germ semis meats Ge pen eee ee SE ee ees ees GE ey mae Gein ios demos eal ieee city eo eee SE ee Coe, eee Ge 


avs 


2 TANT ef, zt, FR 
) 


AK = Tan" Tle ; (224) 
If the Earth were truly a stationary sphere, these equa- 
tions would be sufficient to describe the satellite's 
motion. Because the Earth is actually spinning under the 
satellite, however, angular rates must be considered, so 
that the Earth's rotational speed (shown aS (4%) in Figure 
2.2) may likewise be considered. Toerefore equations 2.3 
and 2.4 above must be differentiated with respect to time as 


follows (note that 1 is time-independent) : 


i) 


DS 


° 


~ ese, 


2 


¥ 2 WP Con Oneasne ae (25) 
Z 
= | x O(z - Z|) LE rT ¥| 


|» 


CoS SINUS | 





Ko, Se (2-6) 

poe @ : 
Substitute equations 2.5 into equations 2.6, 

° 7 CES mGe mn 

sa ) (21s) 
\ Commer SiN od Cent 
Goce d 
X = —— ae nn ee 


. x: 
coS*d + Sin COS Le 
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G@—@sader now that. the scarth is. rotating in an easterly 
directicn; the rate at which east longitude is traversed by 
the satellite is now decreased by the Earth's rotational 


Speed, Or 


; COS Ld 





= Cee 


,) 


\ Ces as ee COS’ 
$ (2.8) 


Ehpioying equations 2.7 and 2.8 for the generation of lati- 
tude and longitude respectively led to the computer- 
Genemeted ground tracks shown in Figures 2.3 and 2.4% fer two 
Meercal Orbital periods (note that the inclination angle is 


60 degrees in both cases). 


7$.0 


. e 
O6 BO - OOF 095 CFE SS CFOS OED OSD OF 
° . 


0.0 25.0 $0.0 


NCRTH LATITUDE 


"23.0 


e e « 
Cee O08 C008 ee SOE soe Oe aE One EEE Oe oe Od OUes emt net 
® 
° 


if : 
« 
: 


e 
t 
3 





Ve 


-X%.0 3.0 %.0 0 &G.00iCiW. it OCIS.) 186.0 «210.0 240.0 270.0 «(3.0 Ck. .0 
ERST LONGITLOE 


Fugmure 2.3 Satellite Ground Track 
Pomerat Hominel PerioG@ sor Ninety Minutes 


73.0 


29.0 


RORTH LAT I TUOE 
0.0 
oo coen ce ceccn ences ofoeee 
| 


eS eweSOeg ert aener a ete eeetO eget edectos cose r ees ee See seSere -PoeeSSGoteesegesesse 


ae) 


: 


Pp 
1 
! 
! 
| 
§ 
; 
i 


O* 26 fee 68 Geet COUT 





* i550 


9c.0 120.0 183.0 190.0 212.0 
ERST LONGI TUCE 


Figure 2.4 Satellite Ground Track 
For a Nominal Period of Twelve Hours 


The iast factor that must be considered to achieve abso- 
iute accuracy in predicting the® sSatelia te*s jorbital Sear ee 
the non-spherical shape of the Farth; to wit, the Earth is 
actually a slightly oblate eilipsoid, which causes a slow 
diash tee or precession, of the satellite's orbital flane 
around the Earth with respect te the fixed stars. pReiss> | 
inciudes a detailed discussion cf the effects of precession 
on the orbital plane, aineluding the serivation, ofeaeto maimed 
For precession rate, which is offered here without DEOCm 


OD 8a. Carew lar JO oie 


° SSG Ne ae L ° 
seis i CoS \e [roa ) 


Mme woe depearss tO be a potentially large number, but as 
Figure 2.5 indicates, is actuaily only a few degrees per 
day. When one considers that the typical tracking satellite 
makes a dozen oor more revoluticns of the Earth daily, the 
@eamwoe 1 the ground tracks of Figures 2.3 and 2.4 would be 
quite small for each revolutior; ‘further, as long as the 
satellite's ground station is aware of the precession and 
makes compensation for it, there is no discernable error 


introduced into the target tracking problen. 
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Cc. EAKTH CCVERAGE 


The key guantities to consider when discussing Earth 


coverage for either a single satellite, or an array of 
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satellites, are the related concepts of swath width and 
revisit time interval. Swath width relates directly to that 
area viewed by the satellite as it passes overnead; revisit 
time is the meaSure of how frecuently a particular point on 


the ground is in view from a passing satellite. 


aL PE Sam ee a 





——— oe i 
[ay ee ee ee) 


Figure 2.6 Geometric Swath Width 


The first aspect of swath width to be considered is the 
maximum viewing angle, or Geometric Swath Width (GSW). As 
Showh in Figure 2.6, this is the maximum viewable area by 


the satellite, and is related te the altitude as follows: 


© Vo 
Cos > = 
hase van ; 


OL 


at ‘oc 
©- ~ees || (2.9) 


Or, in terms of distance, 





ek % 
GSwW= 20, CS as (2.10) 


Some typical values fcr GSW are shown in Table ITI. 
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TABLE II 
GSW as a Function of Altitude 


h So GSW 
150 un cee 1995 nav 
600 na Boia oe 3796 nm 
2250 nm 105. 6° 6338 nom 
5600 nn No Seo 8118 nom 
10900 no 2.3. 9130 nn 
19322 9nn 162.7 9760 nm 
ce 180” 10800 nn 


SS ee ee eet 


ll 


| 


The second aspect of Swath width iS minimum Swath width 
(MSH) recuired for coverage, asa function of the number of 
revoluticns per day. Because tne swaths cut by a sateliite 
in a circular orbit are most widely spread as it crosses the 
Seuator, his simply means ccmputing a, the number of 


eGuator crossings per day, and then swath angle is defined 


as 
— 
oe c V\ (2. 11) 
| 3 
(24 22s)GM 344 cLee Hm ® 
where A: Se (ene 
sO 
5 Sty 
: a Goan 
©, > (1.05 ALO yf | ea) 


Taking into account both the ascending and descending 
trajectories, the typical sateliite orbit can actuaily be 
considered to cross the equator twice per revolution (see 
Bugures 2.2 and 2.4), and so a. is a sufficient swath 
width to ensure coverage along the entire equator, which is 


of course a worst-case approach for achieving coverage of 


a 


the entire globe. This leads to the information summarized 
in Table III; the last (geosynchronous) case is obviously 
unrealistic for tracking applications, as WSW iS greater 
than GSW. 


TABLE IIL 
MSW as a Function of Altitude 


ee | 


| 
1 
| h n Ce MSW 
150 na 16 11.25° 675 nm 
| 600 no 13.41 tae 4 7 - 808 nw 
22 50an i 8 Lo 1350 Sn 
5600 nm 4 Se 2700 nm 
| 10900 nm 2 20° S200 nM 
| 19322 nm 1 180 10800 nna | 
ld 


The overlap, or difference between Tables II and III, 
determines the maximum number of times a single satellite 
can possibly view a stationary target on the equator during 
successive orbits. This numker N can be calculated by 


dividing twice eyuation 2.10 by equation 2.12, or 


=\ Ve ' 
NI _ a oo | | sou 
S Coe 


coe : 





Veo mie = 


Unfortunately, average time between revisitS cannot be 
characterized as 24 hrs/N, as successive swatus will overlap 
on consecutive orbits, and then the target will not be illu- 
Minated for several orbits. At this point it is necessary 
to place additional sateliites in the oreetat pateemn, ween 
constellation, to ensure coverage during every orbit, i.e. a 
Minimum of 360/GSW (xrcunded up) iS required. This will then 


provide n visits per day to the eyguatorial target, where n 
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femeagdain aS defined in equation 2.11. Once this has been 
achieved, insertion of additional satellites wili only 
continue to increase the revisit frequency while correspond- 
ingly decreasing the revisit time. 

The last factor in this tLevisit time problem is’ the 
consideration of horizon-to-norizon coverage along the fath, 
that is, until now the discussion has only dealt with meas- 
uring the time between Successive CPA's (Closest Point of 
Approach) of the sateilite to the target. AS Figure 2.7 
indicates, however, the target 1s actually in view icng 
before the sateilite passes overhead, and remains in view 


until the satellite passes over the far horizon. 
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Figure 2.7 Target Time-in-View for a Single Pass 


If tne target lies directly on the sateiiite's path, 


then the maximum viewing angle © is defined in Figure 2.38 


as 3 \ Gee | Xe 
Sk < a) Cos Tr las + Cas wr \y e 
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Figure 2.8 Viewing Angle 


This translates to a viewing time of 


3lz | ve ee 
ee a Se y i (q+h) ee a ese ie 
*~ S oe ~ : nr . 
3&0 TaGM Lac 


For a surface target and a sateilite at 150 nm elevation, 
this translates to a time-in-view of 8.3 minutes. This 
viewing time will increase as satellite altitude increases, 
and decrease aS the target moves toward the outside of the 
Swath. Regardless, the revisit time between taryjet observa- 
tions 1S always decreased by this viewing time. 

The final facet of ground coverage to be considered is 
coverage at the poles. Havirg used the Equator as the worst 
case for Swath spreading, one might then expect the swaths 


to converge at the poles. Actually, the convergence of the 


Zo 


ground track occurs about the latitude lines of the angle of 
inclination because they are the circles of smallest circun- 
ference to be viewed from directly above; hopefully the 
swath width of the satellite is then great enough to overlap 
the roles. TiS OmaInonm bind! anclination of 60 deyrees, 
the most closely grouped locus of points observed f£cronm 
@irectly overhead are the latitude circles 60° WN and 60° s. 
To see 1£ the swath overlap will at least cover the poles, 
Figure 2.9 must be considered. =A is the half-angle fron 


Sqmation 2.9, so 
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Figure 2.9 Minipum Inclination for Polar Coverage 
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TABLE IV 
Minimum Inclination for Polar Coverage 





h i 
150 no 3.4 & 
600 nm bd . 4 = 

2250 nn SM ew ed 
5600 na py 4.2 
0900 na eo ed 
9322 nu 3 a ae 


> Sa ae 
an |) 


To summarize thus far, a satellite placed at an altitude 
of 600 nm would be useful for tracking purposes; its GSW of 
63.4° and MSW of 13.47 ensure that a target would be viewed 
on at least four (and probably five) consecutive passes, and 
its crbital period of 107.4 minutes indicates that the 
target would be viewed freguently enough to acguire tinely 
and Significant information about its Location, course, Sane 
speed. Inclining it a Minimum of 58.4 degrees would ensure 
coverage of both pcles, and its precession rate of -3 
degrees per day (or akout -0.2 degrees per orbit) would not 
deter its ability to Maintain ~ track on the target. 
Reaiisticalily, a single satellite would not be used; the 
placement of a minimum of six in the consteiiation would 
ensure that no target would go more than one period without 
being okserved, and the placement of additional satellites 
would only then increase the fre,uency of target 


observation. 


D. SATELLITE ORIENTATION IN SPACE 


Having developed a working rodel for the orbiting satel- 
lite, the next phaseset the protlemeyall conmerder ste lcaeer— 


iite's observation of its target, either by passive (tee. 


S10 


TR) or active (i.e. search radar) means. First, however, a 
coordinate frame of reference for that observation is 
required. Thus far the assumption that the satellite can be 
represented by a point mass has sufficed; now, however, it 
must te realized that the satellite has dimension, is irreg- 
ularly shaped, and yuite possiblv undergoes radical attitude 
changes throughout the course of its orbit. It is therefore 
imperative that the satellite be able to orient itselé with’ 
reference to its grcund point, that is, the point on tne 
Earth's surface directly between the satellite and the 
Parcth's center. Thais is because, realistically, the satel- 
Peeomwil! be Measuring (at best) a line-of-Sight range anda 
series of angles (prcbkably altitude and azimuth) to the 
target: it is these angles that must be oriented with 
respect to the reference cirections, in order that the data 
mea to the tracking computer be consistent. 

In considering Figure 2.10, Figures 2.19 {a) . and 2.10 
(e) show the space reference orientation and the satellite 
mergmedsinethat Orientation; Figures 2.10 {b), {c), and (d)} 
Show the initiai position and intermediate rotations of the 
Satellite axes required to achieve that orientation. ite 
should Le emphasized that these rotations are computational 
rather that mechanical; as long as the target falls within 
the sensor's coverage area, there is no need to actually 
Physically rotate the satellite into any new positicn. 
Rather the target coordinates are transformed mathematicailiy 


as follows: 
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Pigure e2-a) Translation of the Satellite 
To the Reference Orientation 


Steps. Rotate around the y axis by il to get the 
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axis ion the "D" direction. 
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Mts tO Got Sromex-v—-z coordinates to N-E-D coordinates, the 


tmensSiemmation is 
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Figure 2.11 Angles Transformed to Rectangular Coordinates 


There is actually another step involved both before and 
Seber the transrenmation of equation 2.13, Hote «OL trans 
forming from the angles to the rectangular coordinates and 
Dac kK. pee aeicldtlcr scrOnR aueinitial Set Of raw anguzar 
measurements to the satellite's x-y-z rectanyuilar coordi- 
hates 1S snown in Figure 2.11. The corresponding eguations 
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transformation from the 
rectangular coordinates to the 
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Meanwhile, the N-cp-D reference 
observed angular coordinates 
is depicted in Figure 2.12, and described matiematically Ly 


Oe i le a 


oo = TANT EL 
B= cos Pla, 
(22 15) 
It is worthy of note that R remains unchanged by ail the 
transformations; only x& and 4 cnange. Thus the complete 
transformation from angles-referenced-to-satellite to 
angles-referenced-to-ground is 
= 
N | | | R Sin (5; SS me 
en cc oF 
= -T \ \ i aan 


te \c Cos Toh 
| (2eaien 


Benet Q = Nr = Wa 
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Pegure 2.12 Rectangular Coordinates Transformed to Angles 
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Figure 2.13 Observing the Target trom Space 


E. SATELLITE OBSERVATION OF A KEAR-SURFACE TARGET 


Eguations 2.16 detail the realignment of raw target 


information along the space reference axes, N-E-D; the 
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Figure 2.14 Satellite-Target-Earth Center Triangle 


final step in calculating target position and velocity on 
the Earth's surface is to translate tke range and the space 
angles into earth's surface coordinates, il.e, latitude and 
longitude. Figure 2.13 shows the pertinent variabies, while 
Figure 2.14 shows the triangle formed by the satellite, the 
target, and the earth's center. The iatter figure presumes 
a surface target; the case of an air target will be 
discussed momentarily. 

To find the earth's surface distance along a great 
circle passiny through the grcund point” “and the Weare jee 
position, one need only find the earth-central ate with 


the law of cosines, as shown in Figure 2.14 , 


oan Ce a ieee ant 
ye fos Car yy = U2 Tide 
Ee CEe A) 


and then the arc lenyth, or surface distance 
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Figure 2.15 Same Triangle for an Air Target 


Note that the angle 2 was not required in this calculation; 
mers Dequired, however, for calculations involving an air 
target of unknown altitude, as shown in Figure 2.15 Solving 


for the target height requires applying tne law of cosines 


again: 
ites 
— | ° = re a ore 
= (<, R) ii DR (Kr hs) Cod cy (2. 18) 
and ey follows from the law of Sines: 
Rsv 
ee ea 19) 


and once again, lone Soy 


In either case, cnce the great circle distance from the 


Seeumd pOlnt tO the target is khown, the target's latitude 
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Figure 2.16 Surface Triangle for Target Position 


and longitude may be calculated by using the spherical 
triangle shown in Figure 2.16 . Using the law of cosines 
for spherical triangles, and realiziny now that it is actu- 
ally more convenient to work with earth-central angies 


rather than yreat circle distances, then 


COS (<tc AL = cos (A6-Ls) Coss T SiN Ctr SiN Cus XK 


Or, eguivalently, 
Ly = Sa Fs Ls Ne t Costs Se Ces rar (2 5210)) 
Sam ileal cage, 
CUS Y= Cos (ae 7) cos (40 ty) + Se OGe ae -\) COSIZL. ae 1h) 


Cl. SOlLVinGeeors jar 


4 | COSD— Vin 7 Sea. 
JC =F Ces ie. = 
Cos le Ces Ls 
(22522) 
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It snouid be emphasized that Ecth latitudes are referenced 
North here, and both longitudes East. TO Switch to south 


latitude and/or west longitude, one need only know 


DOmm@Elat <alu,7meoeLat = —N Lat. 
LOR TOM VU, te roncGg = =F Long ; 
POmea longue, Wo LCG = 360 = £ Long. 


Sample calculaticns usSing these eguations reveal that it 
takes a Significant amount of target elevation to generate 
an appreciable change we mca ederange Of 1000 nm £roa 
feo meellate dt elevation 600nm, a target flving at 25000 £t 
would have af yo angle only (eeemidiaterenteerom tat of a 
Surface target. Therefore, depending on its sensitivity, 


the satellite sensor may or may not distinguish bétwen 


target types, 1f dependent cra singie look at position 
only. The necessarily greater speed of the air target, 


however, will ultimately become obvious if the satellite 
system is either carable of calculating range and bearing 
rates based on multiple observations, or else rate infcrnma- 
tion 1S ah aintrinsic part cf tne single observation. 
Therefozce, this model must frovide for the satellite's 
ability to observe R, oe paid 8 in addition or aS apposed 
SOR, xX , aid a The transformation eyuations 2.14 through 
2.23 must then be differentiated with respect to time to cet 
the range and bearing rates referenced to WN-E-D, and then 
the target's latitude, longitude, and altitude rates. 
Differentiating equations 2.14, 


od 


ae VSN ACsik, + RA CES [COs ume Siw B Shee 
mt eee San Sea. * Q fA, cos B Sin dh, 4 AL ssf Cor hy 
aa c 


3: Aces £ - RB aN 


(2. 24) 


Bis, 


This gives rate of change of the target's coordinates 
referenced to the satellite's axes and based on the target's 
motion; if the satellite can instantly fix its position in 
Space, then x-y-Z Can be transtormed to N-E-D Dy ape ying 
the transformation eguation 2.13 to correct for satellite 
attitude. If the satellite's own rate of rotation around 
its three axes is taken into account, however--1.e. tne 
Satellite 1S tumbiing fast encugh to arfect significantly 
the observed angular rates--then that transformatron ~2 318 


must also pe differentiated, i.e. 


Nw] fe espns 
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In the next step, equations 2.15 are ditferentiated to 
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Having transformed tke raw rate data to the reference coor- 
Sm@Gate ates with equations 2:24 through 2.26, these rates 
acre then translated to earth's surface coordinates by 
Peer etrentiating equations 2.17 through 2.253 . 


Then, from 2.17 (surface target), 








o Zee . 
“ “Ct hs) Se (Ze 
meom 2.18 , 
C R ( R- (ths) ces ) ROG: hs) 4s . 
—. <<. °° °°» | (2. 28) 
Protez. ts (air target), 
¢ («,+lay) Iasi + eB cos B|- Rhy sin A 
as pcm } (Z. 29) 


(Ge no) cosy 


and in either case, Racy. 

mEom 2.20 , 

hi. = _) E Jousts Co5o = SiN saypeesa|- a costs Sif Sus 

tp | ces locos co$L—- Sin i, || | (2. 30) 
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These equations 2.24 through 2.31 completely define the 
derivation of target latitude and longitude rates tkased on 
raw angular and range rate data available at the sateilite. 


If applied seguentiaily by a computer, the target's rate o£ 


travel could easily ke calculated. 


4 4 


This completes the descrirtion of the recorded data 
available from the orbiting sateliite and its observaticn of 
a near-surface target. The next chapter is devoted to 
reconstructing the farget™s motion frome enis Vinee r a en 


this will then lead to the observability analysis of Chapter 
IV. 
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Tif. TARGET MOTION 

It 1S aenceforth assumed that the satellite's fositior 
and orientation in space are knecwn, and thus observed target 
data can be translated into ground coordinates in accordance 
with the discussion fresented in the last two sections of 
G@oapter Il. This Leaves the froblem of reconstructing the 
target's motion from the observed data. With an eve toward 
the observability analysis upcoming in Chapter iV, an effort 
will Fe made to employ a linear target model, while the 
observations will be necessarily highiy nonlinear. The 
alternatives would be nonlinear eyuations describing toth 
the target's motion and the okfservation, Om sa DOM near 
plant with a linear okservation achieved by another cocrdi- 
nate transformation, Doth Of whieh would only further 
complicate the problem at this point. They do present 
viable possibilities, however, which might be investigated 
as part cf a search for convenient coordinates conducive to 
appropriate nonlinear observabilitv analysis and tracking 


desiadn. 


A. TARGET MOTION RECONSTRUCTED FROM SATELLITE DATA 


mn Chapter Ii, aigorithms were developed for trans- 
formemgwraw Satellite sensor information into ground coorédi- 
nate data; given range and tearings, the computational 
result is target latitude and longitude. If range and 
bearing rate information is also available, the result is 
then latitude and longitude rates. 

To reconstruct tke target's motion from this information 
teeta vply strazrgnutforward. The first step is to translate 


the latitude and longitude into linear distances, and then 
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Figure 3.1 Earth's Surface Projected onto Cartesian Grid 


compute the velocity components if they are not readily 
available. If the earth's surface 1S projected upward onto 
a Cartesian Grid With 1ts orig rmyeat ee, d ) = (0,0) as in 
Figure 3.1, then a fosition (2,6€,n,) would be descriapedmay 


equations 3.1, with h. unchanged from Chapter If. 


Y= Vi + 


e@=v. hr costly (3.1) 


If rate data 1s available from the sateilite, then the 
nortneriy and easterly velocity components are computed 
Simnoly by differentiating eguations 3.1 with respect to 


time; cbviously, h.. 1S once again aS imme Chapter [iw 


A = 18, eae 


V 


\) 


J) 


N.2 @ = ary jie | (ee2) 


yy 


Oh the other hand, if rate data is not immediately avaii- 
abie, the next position observation is required, and the 
velocities at time k are approximated by 
ae - nde)-nv-s 
VA ce A Cr) ~ ej 


Nl (y= e(= e(k)- elx-\) 
TT } 


Ux) = Meteo Cu) 
\ 


(3.3) 


: ee: ay 
in either case, then velocity N= | Na *\i 


and course ay = TAN ba 
A 


@e She@Wn in Figure 3.2. For an aircraft flying at constant 
omereuae, h-will be zero; for a surface target, both h, and 
ho will be zero. Otherwise, this first-order model allows 
Lor nonzero velocities in the n, e, and h, directions. Bice) 
avoid the necessity tc include any higher order variables in 
the target descripticn, the mcdel will represent acceiera- 
tuome as the controi vector u(t) in all further discussicn of 


the target's performance. 


B. SATELLITE DATA OBSERVED FROM TARGET AOTION 


Thus far, the result of tae discussion in Sections IID, 
IIE, and IITA of this report effectively has been to develop 
an algorithm for describing the motion of a near-earth 
target, given the parameters cf the orbit of the tracking 
Satellite and the orserved target data. No matter what 
method of analysis 1S employed, however, any observability 
study reguires first a target mcdel, and then some series of 


observations based on target kehavior; therefore it now 
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Figure 3.2 Resolution of the Velocity Components 


becomes necessary effectively tc work backwards through the 
ejuations of those three previous sections and develop a 
system of equations that feature the target as the plant and 
the satellite as the observer. It would of course greatiy 
Facilitate the observability analysis if both the plant and 


observer for the system could be linear, and take the fcrn 


x(t) = Ax(t) + Bu(t), 
y (ee =e Cx (eye 
In that case, the observability study of Chapter IV wouid 


become straightforward, and this aspect of this probiem 
would no doubt have been fully explored years ago. 
Unrortunately, as the complicated angular rotations and 
Spherical geometry of Chapter II portend, there is no real- 
istic way to pick variables sc that the problem wili te 
rully linear. It remains then to experiment with either a 
Jinear plant and nonlinear observer, or vice versa; it would 
ultimately be advantageous to avoid if possible the 


completeiy nonlinear system described by 


n(Ce)) weet) 6 (KE) p18) 
gp (ie) SS BUR (Ae) G2) - 
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Now, the discussion of Chapter II essentially derives 
Pees system equations for describing a target's motion in 
terms of space coordinates; however, to have a near-earth 
target described in terms of R, o& , ., and their deriva- 
tives would not necessarily be useful, as it would require 
all the coordinate transformations of Chapter II to be 
performed over again to relay any useful information to the 
ground. In addition, the presumption of a constant (in the 
absence of control input) course, speed, and/or rate of 
climb (air target) would be mathematically impossible, since 
these coordinates would not be available; tnerefore, since 
assuming constant Re sl pee and v5) would still not achieve 
Gomoeanit V., V., and h_, the second derivatives of R, x , 
and £§ wouid be required, i.e. the differentiation of equa- 
tions 2.24 through 2.31. To avoid this unnecessary conmpli- 
cation, the reverse derivation will be performed here, i.e. 
a description of both target motion and space observation in 
terms o£ earth's surface coordinates. The model will be 
derived for the most general case, i.e. all variabies non- 
zero. The constant altitude aircraft or surface target will 


be considered simply to be special cases of this. 
i. Target Model 


If the target is presumed to maintain a constant 
course and speed in the absence of a control input (acceler- 
ation), a simple one-dimensional model might take the forn 
Gir 


% (4 = -: a) 1° ue Lt) , 
Oe: i 


L 


where u(t) can actually be either system noise or an accel- 
eration. Since three mutually perpendicular coordinates are 


reguired to describe the target position, the system should 
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have six state variables, derined by the guantities devel- 


oped in Part A of this chapter as follows: 


X, = N = Northerly distance travetcd sia, 

X,= n = \ Vcos = northerly velocity (kts), 
X3 = e = easterly distance traveled (nn), 
X,=e = Vz = Vsingd = easterly veiocity (kts), 
X,- = hy= aliteudewapove cactii"s Suerac ec mem, 
oe h,= rate OLvelmm: {kes ; 


and negative values for these variables represent southerly 
Gistance, southerly velocity, westerly distance, westerly 
velocity, (no negative altitude), and rate of descent, 
respectively. For a surkace target, x-Sana x are ezerc, oe 
a constant-altitude air target, x 1S zero. The target is 


thus modeled continuously by 


O\ © CO Ge © oO © 
© © CFG. © 1 0 0 
= GG © \ Cle Gaore 
oie SX C&\-+ CON 
(4 Oe Se eo) % Cel Ca © 24, 
QO 90 CO | ego Ss) ee 
(3.4) 


where u(t) represents either noise in the n-e-h.system, or 


T 
else acceleration in the three directions. Realisticaily, 
because computation reguires a finite time and one satellite 
cannot normally be dedicated sclely to observing a Single 


target, the observations will necessarily be discrete, and 


nopefully periodic in time T. The discrete model analogous 
to that of equations seers 
1 1t © Guam Mias) 
OQ \ S Giese TO Oo 
KKHN= lo og i to Of &CK)e loTHRO ub Ck), 
OG © =a CO tm - 
Sc © CO Semi. ome oul 
& Oo ClhaRCm Oko a ‘easy 
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Assuming that external nonlinear forces such as wind or sea 
drag are either negligible and/or unpredictable, this 


completes the target model. 


2. Satellite Cbservation 


TieSobsciveeswiatrix sh Will obviously be a nonlincar 
operator on the state vector. To begin its derivation, the 
taco scoordinates must be transformed to kat 1 tu de= 
longitude-altitude coordinates by solving equations 3.1 in 


reverse , and then differentiating as follows: 


Los iN o 
T t. ees 
eC. © — Ces + e\s Sen 
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Rae eaa GT 
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The ultimate goal of this transformation will be to express 


the srace observations in terms of the state vector, i.e. 


us Q iW, Ox) 
12 R He (0) 
4: mK i OX 
aie = 5 _ 3 ~ — A cs 
= as A Hed} 7 
oe E Hy (x\ 
" B He (%) on 


First, however, it is useful to caiculate in terms of the 


state vector the intermediate variabies 7) and , as they 


were defined in Figure 2.16 . Solutions for 7C are 
teivial: 

T= Ac aAs ) 

He eae. (3.3) 
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To solve for Vis a / , eguation 2.218925 rewrieren sas 


LOLLOWS: 


Cos =Si ly Sit Ue ces (7 Cosiades ge 
(3.9) 


and then — SNP = = ( Cos Lz Samly — Sinsl> COs Cvs 7) 


+Ls (s\nitr COS = Ce sie ies cos) JL cestz Cosky Svs PC , 


¢ \ 
iz a 
: i ~ s1N 2 E (site COS Ly CES 7 = cos Eisai) 
ae (cosq Sinks COS7L = SINL; Co§ Ls) a 2. Case CO sna (3. 10) 
Enough information is ncw available to solive for the 
space coordinates &, (7 5 and their derivatives. ro: 


obtain R, the law of cosines is applied to Figure 2.15 as 


follows: 
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TOMObDtain a, the iaw or sines for spherical triangles is 


applied to Figure 2.16 as follows: 
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To obtain Vow the law of cosines iS again applied to Figure 
yeelo, this time yielding 
re é ra E ie 
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jis obtain J, consider 
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See eh (3. 16) 
Cea (cm) SING 
The ofservation model iS Summarized in Figure 3.3. Moy? 


realistic computational procedure would of course use the 
intermediate variabies L, x Ds ah 4 , and their deriva- 
tives; for the sake of completeness, however, it siould be 
noted that the observer functiors could Le written as func- 


tions cf the original state vector. 
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Figure 3.3 Summary of State and Observer 
Equations Including Intermediate Variables 
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3. Iwo Simpler Cases 


Pre Ssrgnit 1 cant compiexity of these eguations 
Suggests that some Simpler representation of the systen 
might be useful in the observarkility analysis, Di wo nly Ee 
test the method prior to its employment onthe system of 
mgure 3.3 . One aprproximatior might be to presume for the 
moment that the sateliite is directly overhead the target, 
and that relative motion between them iS arbitrarily small 
enough that tne Earth 1S approximately planar beneath it; 


Hone the trigonometric functions, 
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Figure 3.4 Satellite Directly Overhead the Target 





If one considers that fcr a satellite directiy over- 


head the target as sShcwn in Figure 3.4, 
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Figure 3.5 Perturbation About the Observation 
From Directly Overhead 
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Figure 3.6 Resolution of Beta Along the N-E Axes 
then fcr a satellite nearly overkead the target, the probien 


can be considered as a perturbation about the directly over- 


head case. In Figure geo, let 


54 


en cca 


eee = Ni ee 
eee AC = Nee see 
x = Ae = Ne = NC. Ges tk Ks, 


Ms a Al = Ws- nq , 

% 7 At —he . 
To fully linearize the problem , it turns out that 6 nust 
pe eliminated; 


1a ae 


) 


this can be done (in this special case only) 
as shown in Figure 3.6: 
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and the corresponding discrete system plant is 
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with the same observer matrix as above but with values taken 
only at discrete points in time. 

Experimental computations show that for a satellite 
at altitude 600nm and a target at less than 60° latitude, 
these approximations yield abkout one percent error or less 
for vo: less than 5°; however, the maxinun possible over- 
head time during which this model applies 1S about 20 
seconds. It follows then that this sort of observation 
could only be made a very few times during the satellite's 
revoiution, severely limiting the real-world usefulness of 
this linear model. For the purpose of examining the observ- 
ability of the target under these conditions, however, the 
modei should still prove guite vwseful. 

The other Simplified model that might be instru- 
mental in gaining some insight into the method of the 
observability analysis is the simplest nonlinear case, i.e. 
the satellite tracking a fixec foint on the earth's suriace. 
Simplifying the equations 3.5 through 3.16 yields the system 
depicted by the equations in Ficure 3.7 


This completes the various derivations for all the 


target models. 
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Figure 3.] Eguations for Tracking a Fixed Target 
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This Chapter is subdivided into three major sections, 
specifically the okservability analyses of the three 
different target models develored in Chapter III. Pirst, 
the linearized model will be analyzed by the method of 
[Ref. 2], and then the stationary target and the general- 
case, nonlinear, moving target will be analyzed by the 


method developed in j[ Ref. 1]. 


A. OBSERVABILITY OF THE LINEARIZED MODEL 


For convenience' sake, the linearized model has teen 
subdivided into four cases, nanelvy surface and air targets 
with discrete and continuous nodeling for each. The method 
of [kef. 2] is the formation of the observabiia: y pact 


from the system 


x= Ax + Bu 

hes 
such that G Se Fes. Ue ee _— ae | 

~ oe awe ' vw ~ } ere aa 
for n observations, and for sys tenpobeerva bits ys lee 
HouSsil Agu tag. In the following derivations, the mininua 
amount of information necessary for Oobservability is 
considered. 
1. Surface Target Continuous Model 


The system modei is 
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For two or more looks, 
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Pemethe next observation, A‘ = 0, so tnere 1S no need to 
take the process any further. If only position observ- 


ability is desired, rather then full state observability, 


then 


2. surface Target Discrete Model 


Now the system model is 
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looks, i Oo oC = | 
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is reguired, and 
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which sinplifies easily to the same nonsingular matrix 
obtained in the anaiogous continuous case. 

Obviously for this simplified case, the transition 
from continuous to discrete modeling made no difference; it 
has been implicitiy assumed that the discrete samples are 
taken often enough that all maneuvers (i.e. acceierations) 
are observed, 1n accordance with the Sampling (or Shannon 
Information) Theorem familiar from communications applica- 
tions. Otherwise, observability will decrease significantly 


as information about the target is lost. 
3. Air Target Continuous MNcdel 


— ae [=a 


The system 1s now modeled by 
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is required, and 
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EOE complete obServability. AS in the surface case, if only 


position observability is desired, then 
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is adeguate even with only the single observation. 


4. Aix Target Discrete 
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is reguired, and as expected, G = C “again 2) For tues aes 
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which simplifies to the same ronsingular matrix as in the 


the assumption has been made that 
Nyquist) rate, 


continuous case. Neal, 
information is sampled at the minimum (i.e. 
which is twice as often as the target's highest freguency of 


Maneuver. 


Be. NONLINEAR MODEL CF THE STATIONARY TARGET 


Applying the method of [ Ref. 1], summarized here 


briefly, Concerns primarily the analysis of a Jaccbian 


Matrix defined by 
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for each observa tlonsotean n-sy seer. A nonzero determinant 


will show that all the states ir X are connected to at least 


one of the measurements y, and So an inverse function exists 


such tkat 
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a= wl Cy) 
Weeeovseivdbality, however, this mapping of y back into x 
must ke one-to-one, so that a particular series of measure- 
Bees y wt atways detergane a wnhigue initial state x (t,). 


This requirement is satisfied if a second condition is met, 


Lecce Detuste ve POsttivemmer negative definite, and the 
Mapping unigue with finite coverings. The complete procft of 
these concepts is presented in [Ref. 1} and {Ref. 7}. he 


multiple observation variables (i.e. numbering m), ann x mn 
matrix is formed; if any n x n matrix can be extracted that 
Satisfies these conditions, the system wiil be okservaole. 
Por the Six-variabie, six-observation case upcoming in 
Section C of this chapter, this will lead to a great many 
possibilities and problems, but for the two-variable case of 
mais section, the method should prove very usefui. The 
approach will be to first consider each observation sepa- 
rately, and then in combination. 

Starting from the eyuations for tne observaticn of the 
Sedi~ronary target dericted in Figure 3.7, derivation of the 
Jacobian matrix for y = R is accomplished by means of 
MRovicit differentiation along the lines of functional 


dependency depicted in Figure 4.1 . Consider 
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Figure 4.1 Functional Relaticnship of R and R to n and e 


Thus the Jacobian matrix for R is written 
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Simplification of the determinant results in 
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which is obviousiy ncnzero in the general case but mav well 
have muitiple soiutions. Based on the criteria just 
defined, this indicates that bcth states in x are connected 
to tne observation R, but the connection may very well not 
be one-to-one; some further analysis is obviously required. 
MiewletcEMiidilts Oleic J-Matrices for @edch of the other 
observed variables may well prove simiiarly ambiguous; 
considering the observed variables in combination and 
attempting to establish nonsingularity and positive/negative 
definiteness of the resultant wyriad of Jacobians may well 
prove diitficult in the extreme, if only inthe sense of 
requiring a copious amount of ccnfusing algebra. While the 
basic geometry of the problem as developed in Chapter II 
leacdS one intuitiveiy to the conclusion that a Surrace 
position can be uniquely established if either (XS ) or 
(R, % )--but not (R, 8)--data apecmoy ct lapremmdemonstrating it 
by this method proved a nammoth computational undertaking. 
Needless to Say, when the proklem is expanded to the four 
variables of the moving surface target and then the six 
weanpaables of the air target, the algebra worsens several 
times over. The only alternative seemed to be a transforma- 
tion of coordinates, Specifically toa rectangular systen, 
so that the differentiation could include some polynomials 
instead cf only trigonometric functions. The transfcrmation 
is accomplished as fcllows: establish a plane at a distance 
oe tO (ne -h{ js below the sateliiate perpendicular to tne 
vector connecting the satellite to the earth's center, 
determining in erfect the target's horizon. For a surface 
target, the point of intersection is the satellite's ground 
point; for an air target, the intersection is h above the 
ground pcint. Then the target's position 1S projected up 
into the tangent plane, as shown in Figures 4.2 (a) and (b). 
Tt is now possible to eStablisn the relative coordinates as 


follows: 
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Figure 4.2 Projecting Target Position 
Into a Perpendicular Plane 


As every set of true target cocrdinates yields a unique set 
of relative coordinates for a particular satellite, the 
mapping into the Cartesian coordinates 1S one-to-one, and an 
oLbservarcle target in one syste will Fe observable in the 
other. The complete set of state and observer equations in 
this coordinate system 1s as illusteated ingpijures. (Dd) 
and Summarized UnorPmgure vias. For the special case of the 


Stationary surface target, the system simplifies as follows: 


66 


mince Only position 


will ke 


>< 


aie 
M. > 


1, 


oN C8) 
ae Ls 
eC Adm i.) Costs 
aN Ag coslky a L Neg ( ee \\ SiWhk, 
hs 
O 


-~ 


The plant equation is unchanged, and the observations are 
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Figure 4.3 State and Observer Equations 
For the Cartesian Version of the Problen 
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which 1S again nonzero in the general case, and which will, 
on refiection, lead tO a Whigque Solon Or anys eaeewcr os 
target and a specified satellite; in fact, equation 4.1 


Should then also yield unigue sclutvons, Dee just eae 
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Sie = O-One COErrespondence Gfethe Coordinates and the observ- 


ability froblems, but also due to the physicai reasoning 


ierasthated in Figure 4.4 . 
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Figure 4.4 Tracking the Stationary Target 


The sateilite's ground track is a curve, and there is 
Only one possible foint of intersection of even three 
circles whose centers lie upon it. PittseOhesemvant tity wild 
Ee more clearly demonstrated if tne single observed variakle 
to the stationary target is ~ ; obviously, the intersecting 
lines oi bearing must define a single point. Further, since 
for a surface target every R defines a unigue 1S) and vice 


’ 


versa, the A-only case should ke observable also. 
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Consider first the x%-only case: 
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The numerator is identical to that of the J, determi nda 


Figure 4.5 illustrates the situation here. 
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Figure 4.5 Azimuthal Bearing-Only Track 
Of the Stationary Target 


For the record, in the S-only case, 
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So again the same fcrm of determinant will be achieved, 
eeeeit wath a different multiplicative constant. Thus for a 
target known to be stationary, the observation o£ any of the 
@adntities R, & , or yG is enough to fix its position. 
Presumavly the same result couid have been achieved in the 
Original coordinate system had the algebra been carried 


PrEeough . 


C. NCNLINEAR MODEL OF THE GENEBAL TARGET 


Once the target is permitted to move, or course, the 
unigue-intersection arguments just developed no longer 
apply. The eguations shown in Figure 4.3 now are pertinent, 


Piteed is defined by 
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The derivation of the Jacobian matrices and their deterni- 
NantS iS now a Major froject ir itself, but it follows the 
Same pattern previously established. In the interests of 
brevity, the requisite algebra is exciuded from this discus- 
Sion, but included in this report as Appendix A. After 
Smet icatven by =rew reduction, the J matrices are as shown 
in Figure 4.6, with the ) and 4 terms as specified in the 


appendix. Further Simplification of the determinants may 


71 


well be possible, but these representations should suffice 
here. It 1S obvious from the occurence of the zero rows and 
columns that all the determinants excluding those of J - and 
J, are zero; computation of these last two determinants 
proves them to be Zero also. Considering tne observed vari- 
ables in combination, it would seem that some variable in 
addition to R or R is aiways necessary for observability, as 
they have the only Jacobian matrices with zero rows. YLhile 
the possible combinations of six columns out of the thirty- 
six are so numerous (approximately two million) that it is 
virtually impossible to explore them all, only one combina- 
tion with a nonzero determinant and positive or negative 
definiteness is required for observability. 

What one would expect, based on the insight gained in 
Chapters II and III from the mathematics of the derivation 
of the model, is that any combination including or e and 
at least one poSition variable (i.e. R and ee ao and oe 
etc.) should be observable, and that any determinant formed 
from a combination of columrs from the appropriate J 
Matrices should prove it (excepting of course the zero 
columns in the «& and od matrices). The combinations of 
Pee i a B shculd however prove unobservable for the 
reason of a lack of positive/negative definiteness, that is, 
an ambiguity of target track resulting when reconstruction 
of the initial state is attempted from the measurements. To 
prove this Condition true for each of the approximately 
130,000 possible matrices 1s however obviously anh unreal- 
istic undertaking. 

One Slightly simpler case cf particular interest night 
be that of the surface target, which involves the first four 
rows by four columns in each matrix. If the accelerations 
are nenzero, all the Jacobians except R have nonzero deter- 
Minants. The el ae, and 8 determinants are not posi- 


tive or negative definite, indicating that there is still 
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Figure 4.6 Jacobian Observer Matrices 


For the Six State Variables 
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some ambiguity in the observaticns and so the target is not 
@bservaciec. In the case of the o Jacobian, showeve: eee 


positive definite, i.e. the < terms drop out and 
\3, = (ea. -nade) (4.5) 


This indicates that a maneuverirg surface target is observ- 
able when azimuthal t$earings only are observed; this case 
corresponds very closely to the two-dimensional bearing-only 
tracking problem long familiar to submariners. If the target 
does nct maneuver, then the Jacobian becomes aon- 
Singular and the target is unofservable, and then more than 
one variable of observation, to include X or ol ; is 
required for observability. 

It becomes obvious why the observability of nonlinear 
systems has long been considered a highly difficult problem. 
At least as much insight was gained into the probiea from 
its mathematical derivation as from the application of this 
Method of observability analysis. One area in which the 
method might also be useful 1S in evaluating the sensitivity 
of the svstem observability to fluctuations. Specifically, 
the terms in the J Jacobian are each smaller than the xX or 
R terms by a factor of 1/\nt+e’, which would, in the general 
case, mean a difference of three orders of maynitude. T has 
would certainly lead cne to presume that the system would be 
more sensitive to fluctuations in 4 than in R or x. 

To summarize thus far, the geometric derivation of the 
coordinates in Chapter II leads one to believe that ail 
three coordinates R, &% , and fb, or the derivatives of no 
more than two of them, are necessary to reconstruct the 
taryet motion in the most general (air target) case, wnile 
either the (R,X ) or (XH, 4) combination is necessary and 
sufficient for the surface target. EXamination of the 


linearized model in this chapter tended to validate those 
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insights; then the examination of the stationary surtace 
target led to the coordinate transformation, and to’ the 
conclusion that any position variable of observation would 
surfice LOF observability of this type target. 
Consideration of the yeneral target model showed it to be 
unobservable for any Single variable of observation, and 
then a brief examination of the movinj surface target led to 
the important conclusion that ar c{-only track is possible, 
at least if the target maneuvers. Further, it appears that 
the system is far more sensitive to fluctuations (i.e. 
noise) inthe 4 measurements than R or Q. In Chapter V 
some cf these perceptions will Te tested by the designing of 
a deterministic nonlinear observer and then tracking tne 


target with an EKF. 
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Ve. FILTERING AND TRACKING 
The objective Of@this Ghapter 1S tO dpply the veamiicuc. 
Chapter IV. Specifically, this will entail deSigning a 
deterministic nonlinear observer, and then an EKF, in accor- 
dance with the methcds discussed in [Ref. 3] and [Ref. 8]. 
The original coordinates of the problem as developed in 
Chapters II and JII will be employed, rather than the 


rectangular coordinates of Chapter IV. 


Ae. DETERMINISTIC NONLINEAR OBSERVER 


Actually, in [Ref. 8], Kirk only discusses the design of 
linear okservers, but by emplcying the now-familiar tech- 
nigue of linearizing about an eStimate with a first-order 
Taylor Series expansion, the technique can readily be 
extended to a nonlinear system. Some alternative methods of 
nonlinear observer design are discussed in [Ret Sap 
[kef. 10], [rReft.. 11], and [| RetimiZie 


For tne system under consideration, 
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Assuming an observer-estimated trajectory x, the observation 
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Figure 5.1 Plant and Observer for the General Systen 


Ignoring the higher order terms and assuming for the mcment 
a constant velocity target, the linearized system can new be 


descrited by 


(3-3) 
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Now, Since the methed of [Ref. 8] assumes that only the 
measurements y are available, the plant matrix must be simu- 


lated by the observer to generate successive estimates, i.e. 


a RS (ee } 


\>> 


However, unless the estimates are coincidentally initialized 
to values identical to tue initial conditions of the actual 
systen, there will be an inherent propogation error. 
Therefore a correction term must be added to equation 5.5 as 


foilows: 
A rn rm 
x = Ax GH (x) —- fipey (5.6) 


where G is the observer gain; the @deteruimation (of@Giwa ta 
discussed snortly. Faeuations 5.7 through 5.6 Can be J1igeS 
trated by the block diagram of Figure 5.1. 

To evaluate G, the first step is to substitute equation 


5.4 into equation 5.6 so that 


x? « 


a 5 ae 
eae ol (Ox) ; (5-7) 
y ox ‘ — = 
It becomes useful to define the error in x as 


A 
A ck 


ie 


) (5.8) 


54 1s defined by subtracting egtation 5.7 from equation 5.3, 


Lawes 


Substituting eguation 5.8 where applicable, 
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a or 
Or = Fx Where fF = A=G A\. Pa S29) 


ine solution to equation 5.9 is 


Et 
G(t)- e 7 X(c) (5.10) 


foremecetverges fOr any sb with regative real eigenvalues. 
fee me oat ay detainee ine ak Case Where H 1s a constant 
Battig, a unique G Matrix exists for any desired set of 
fegeeave cigenvalues or Ff; however, for the general case, 
Sih e |, obviously varies with state and time, G must 
generally be a time-varying or state-varying function. Thus 
cw EOulem remains, how to calculate G. Many methods have 
been demonstrated for the case of the time-varying linear 
observer which could be extended to cover this linearized 
case; the method of choice here, however, will be the KF, 
Wilch reguires some ovreliminary discussion prior to its 


employment. 


B. EXTENDED KALMAN FILTER-- BACKGROUND 


Poche. and (ref. 8] discuss the analysis origi- 
nally employed by Kaiman and Bucy to develop the linear case 
and the so-called Extended Kalman (nonlinear case) £4.) beim 
It suffices here to state that the theory applies to 
stochastic systems, rather than the perfectly deterministic 
system considered thus far, and invoives perturbation about 
ees iZ-dgeestimate Of Ehe target Erack in the same 


fashion as Section A above. 
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Figure 5.2 Plant and Observer for the Stochastic Systen 


Figure 5.2 illustrates the addition of measurement and 
plant (or "system") noise to the probien. The describing 


equations of tne plart are now 


Se ae (5.11) 
dacs ~~ oa ae | 
| eS) (5.12) 
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mete W ts "system noise, and vy 1S measurement noise. 1B 
important to note that the form of the observer nas not 
changed since the discussion in Section A of this chapter; 
now, however, the filter analysis provides a _ methcd for 
Geecueldting Lhe Sequence Of gait matrices G. 

That last statement concerning diese dGedtees Of —=G"Ss 
Suggests the convenience of a discrete (as oprfosed to 
continuous) representation of the systen. COWstaeri nc Ghat 
the iterative calculations required to compute the succes- 
Sive estimates x will require a measurabiy finite time- 
increment, it makes sense to sork with the discrete systen 


equivalent to the continuous system discussed thus far, 1.ée. 
x(x+1) = Ox(k) + Lw(k), Gaus! 


¥(k) = H(x(K)) + ¥ Ck), (5. 14) 


_ 
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and w Cx) = pe We (K\ 2 N(c,5,7) , Sic =e = a =o(a5 


We CK) 
Wah) 
Ck) = ck) , ui Ox] Pe Al(¢,5,1) sn Ti - ae =O). 
Vi Ck) 
Vg Ch) 


That is, both the system and observer noises are presumed to 


be white independently normaliy distributed with mean zero 
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and specified variances, which are independent of the 
Gaussiar random’ variable x(0)- The observer eguation 5.6 


Wiil take the forn 


¢ r 


UES) = q aCK) + GC | AOxcel - H (200) Gass) 
It is now useful to introduce the double-subscripted nota- 
tion required to compute the G(k) matrices in the EKF algo- 


Een, 1L.es 
A 
X(KjkK-1) , 


Which signifies an eStimate at time increment k tEaséd on 


information available at time ircrement k-1, as opposed to 
“A 
x(K 1K) 


which represents the estimate at time k given the informa- 
tion available at time k. 

The final element reguired to understand the EKF algo- 
rithm is the concept of the elror covderancamy mace Cuream 
whick can also bLe thought of as the expected value of the 


error vector squared, i.e. 
P(k) = E(()"). 


Its derivation and significarce is fully explained in 
[jRef. 8], while [Ref. 3] discusses at length its initializa- 
tion. Suffice it here to state that it iS an intrinsic part 
of the EKF algorithm, which can now be summarized. 

Given the state estimate and covariance matrix at time 
k-1 and the measurement at time k, the state and covariance 


at time k are estimated as follcws: 


(6 1s) 


J 


fckik-) = & (er[er) 
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P(kIK-1) = @ P(x-r(u-r) 4 al al, (5.17) 
where (CUES Es SS 
iz eS <x G 
ome Gr 


the EKF gain matrix is then calculated as foliows: 


Z| 
G(k) = i (x \u-) WOW LO P (Ai) (4) + a (5. 18) 


oN | 
O# Malin 


| | Weare: 
woere h (k) = AND {2 =: ae 
Ve 
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oO ND 
Ga © 
Finally, the state vector andthe covariance matrix are 


updated as follows: 


R(kiky = \ (4( | + G Cr] Roe CX Cishe-n)) (5. 19) 
Bash): = ue G Ck) Wc e (x\%-1) (5. 20) 


Mmeiethat £emains 1s the question o£ initialization; [Ref. 3} 
develops fully the method used to initialize both the state 
vector and the covariance matrix, specifically the use of 
the first two target observatiors to generate initial esti- 
mates for both x and FE. This method is illustrated in tne 


specific example that follows iagmediately. 


C. TARGET TRACKING WITH THE EKF 


This section discusses the development of a computer 
projram designed to implement the basic Kalman-Bucy filter 
algorithm discussed in Section B above. While admittedly 


Geit@iem Iteserves aS a first-cut tracking algorithm and can 


om 


be used in conjunction with the observability analysis of 
Chapter IV. The proklem divides itself naturally into two 
parts, namely ceveloping models for the satellite and the 
target; and then using the EKF algorithm, that is, equations 
5.16 through 5.20, to see if the satellite can track the 
tab get. 


1. Satellite and Target Models 


This is the time when the various equations devel- 
oped in Chapter II (satellite) and Chapter III (target) are 
necessary. For the purposes of this problem, a satellite 
with an altitude of 600 hm and an orbital sinclinationmes 
63.4 is used, as this is a tyfical orbital pattern used to 
ensure whole-earth coverage in a multi-satellite constella- 
tone A surface target simplifies the probiem slightly, as 
x need only be four-dimensional; thus the target chosen isa 
surface vessel on a heading ci G75 ake é speed of 28.3 
knots, which equates to a northerly velocity of 20 knots and 
an easterly velocity of -20 knots. At time zero (i.e. k=0), 
the satellite is located overhead the point (0 N,O £), while 
the target is located at (45° N,39 E), which is initially 
beiow tne satellite's horizon; the program is designed to 
compute (R, o& , G ) as (0,0,0) until the target comes in 
view. The computer program listing appears it Appendix B; 
of perhaps greater interest is the data which appears in 
Table V, which shows the program run for 23 minutes with a 
time step of 30 seconds, i.e. kK running efron 0) to 5c.) on 
this single pass, the target was in view for 117 sterfs in 


time, or 19.5 minutes, a fairly representative viewing time. 


2. EKF Operation 


In the main, the EKF algorithm used in the computer 
program that appears as Appendix B 2s the dprect  apprmcarien 


or the equations S2166 through Se20 to Stewart cuemia: 


84 


Rr re ee Oe en 


Vv 


TABLE 


LaaOQaooCoCcoe@e er inet 71OS OTN SO OMDMNOWMT (0NS SOCOM Ver Cot caooEee 

HOD DODO DO DODO DMN FANOMH)ANANDMNANHNOSFNVYO DO AMON FOOMMNDBDDDMNOANMODOOO0O0O 

oo 

NQNODODOVDODODVDODO DODO MM ON ANKDOMNDIFOFrAWDANDN FOTMDIDONMHNVUOFE DWDDDOOCOO 
MAW LOD WIN WN D LAY) FS PNA HN MO AT AT TNWMMN MW) 


LeOOoOooaOoOO COON Men COOMNE NMOMNOs OW FROMM OWMKANONas SNC cCcoCceo 
OOo SVOOODOOO OG DONO OOM Ss Ss WN NDOT KH INOMANS ANSI NTOONOOoNyONeoeee 
oF ee ee ee ee ee ee ee 2 i 2 2 2 | 
FADO COD DQDODVO DA GNHOO TANM TNO AHDMNDs NO COP DMPomMoWNstrosromr Coocoo 
ao AMMMMNAMNMN NANNY ttt ATOMS TUNDOSHKANANNANMMM t+ tTNOOWN 

Tm MANANANANNAANANANANANAANSG 


BODDODODOCDD DDO ONN st DWNNOO FSM THNMMH HK KVDODHNDWNOM MONDO ST TrHKOMKRMNWOM OOOOCO 
oe ee ee | 
ZDOOOOO CO OC SOO ONINO mins CEM DNWNWOt Ss COCO ONMK—SOONN Si sNoonooeoego 
cL DADOOOM ANAM FUE DAQNANNA TeHMME ADO SNTONDOMME MM Ww 
ow DONO OM WOlLN AMA HONDO OWOVOMOUW OM PONDOSKANMANWCE ODO 


COON OVO CC COP FFF OVO OWMWMW WY SF FAM MNMOAANANANANT © HOO ODOMDAAOG CO COP FF FP OW WO 
™OODNNANAADAINAADAAHAAHHAAAAADAAAANAAA HAND OAC AAG OOVOV ON OD ON CO 0D CO 00 CO. 00 00 CO 00 CO CO COCO CO CO 
ODOeteweeeeeeee fee © ¢€ © @ © © © © © © © © © © © © © ew ew ee ee ee ee leh hl ehlhlUh hl hmhlUhr hl hmhlc hlUh hl hl hl hl hl lh ll 
AD DANHANANAHAANAAAAA ADDN ANAAKAAHAAAAAHAAAAAAAWA AA AN DVN OY NON AOA ONOVOY 

MIMINAINON ON QUIN CION CI ON AION CY ON ON ON ON ON ONION CI ON ON ON CUIEN ENON AION ONION GION ON ON CIN CGI ON OION ONION OI ONION CON ONION ONS 


FADD PET ONANNANNANNMMOMN St Sat TWN OOOO IEEE DODDNDDANVNNWDO DTS TENANANANMM MS FF st TINO 
(Off Ee 1 He XG | > Meas IDM cep] ee ¥ a> J ae Yaa Yee I ap <> | a> | a cas Yam Tau Y > Y de>) J cee | ae ca Ta Icy | cape a aE aad dk rl eel al el ed ied a Be dh al ne a eas St em el ard 
ro re ee ee ee ee er SS a i OP OO 2 Pe 2 2 | 
LOA WDA LOALAIN LAA LA NA LALA LA WIA LAD LALA LAL LAW ALA LAN ALA LALA AN ALIA LAWL DLA AA LAW A A A A A WI LALA LAL ALN 
SSS EDP af a pt as pms pam me cet ms pon pean fo p rm mm a pms Bam giles eam gies GS et ly re ny gin pe (FOR es a ees em ey Gms nme A eee eS en RM ley gO Get A CD BE gS Je mn me ge 


VDOMMNDNQNOM TON SFMNWODINANKH MEM NMNONNYO TS DOODNHNUONDNNOWNIN DOWNS PONTNNOO SFO NOMEN 
Pe Nae Oia me NT NO ao ON Nee SO Te SO ae 
OQOeeeeerseeee * 6 © © © we © 6 6 6 +e es 8 8 e* © © e @ @ © © © © ee © & © © Be ee eee ef «© e @ © « 8 
BOOK RK AMMSTINNURK ONOOKNM SNOOP WOK NMTOR MONS ODONSTRONING§ Ft ORTONION 
I TO (NN ANAQNAQNQOHMONMN OST SST STAM MNOWOORE OOOO) 


HOD ODNNDOOMOAMNTNOMNOOTVOTSRHDDOO FON SATANDANDOMUONODM DOOMNENNOMN=OMNMOW 
LOW OS DANS HANTOMOME NAOT NDMNHDANAMNODMNCO AS FODHVNOTTTMODMNNOMOKHANOOIN 
oe i Pe 
DTmMNTIE DOTMAOHNONANMNN ODDO AMNOM DONAMST NE DNDTMHANOM HH DNS MNANNANANN 
3 TM SH TMANNNANNAN NAMM MMO MM & St Tt at FF TNM NINN MUI MN OD OD OW OO OOO 
a har a ata Sk aD 2 al ere calle he carle l g ceaaamts crinSt g pe 
i PZ 
BS A EI OC ON Oe cise He Nae tatoo 
TTT TO TMM TANNA ANGINANQVNAANNANANANG 


MOT AMSNOMOADT AM SING QNOTAM SO DAOT AMINO OAT IM TINOH OHM aman 
THEE TAN NANANNANNANAMMONOMMNYAOYYNMNSASt tt sa TST ST NWI 


| SS OD oy OD ME EY mea A pee A EE TD SS LS 


85 


probleng, complete with the intermediate calculations 
required to evaluate some fairly involved matrix multiplica- 
tions and inversions. However, two aspects of the progran 
require furtner elaboration. First, the filter mnisereli Za 
tion, to which Section B of this chapter alluded; and 
second, the computation of the Jacobian observer matrix h, 


which in this case is anything rut trivial. 
a.) Filter ini tiateezaeen 


The methcd develorfed in { Ref. 3] has been 
followed exactly here, thus a two-point initialization is 
used. Specifically, the northerly and easterly positicns 
calculated from the cbservations (R, X, GB ) at time steps 0 
and 1 (i.e. the first two nonzero observations) are the 
initial estimates for x, and x20) ~and from then@ aie 02a 


oldie ies “ard opin me 

Ke) = (Ra, cle, Blo), 
ee = uF (el, dl, B09), 
ZO) = a (Ray, <0), ACY), 
A (ed), £0, 60), 
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ame 
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1) 
ye 


\3 
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on™ 
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paar A 
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a4 
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and Xo (\\ 
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(ge 
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CG 
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Yay C1) 7 (32.21) 





Trom these calculaticns, it Can §readidy Dewsnovnee nde 


initial estimate of the covariance might be 
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The details of this calculatior, which involve taking the 


expected vaiue of the error vector squared as discussed in 
f Fef. 3]. While it 
other methods of iritial covariance estimation 
this method should 


of this analysis. 


Section B above, appear in could be 
argued that 
are available, suffice for the purposes 
MrtsweComeleotesernie ANitwalizatzon of the filter, 


and from time step "2" (third nonzero obServation) onward, 
meme L21ter LUnS accOrdincg to tne aigorithm of equations 5.16 


mmeough 5.20 . 


be. ~OCbserver Function Matrix 


The observer functicn h is defined for all times 
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Computation of the 


(Ap 9endix B) 


derivation. 


is performed in 


Consider first the 


from Chapters II and III 


matrix elements in the 


computer progran 


accordance with the following 


following set of equations 


that summarizes the computation of 


(Re X, J) from x (surface target case): 


on 
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Figure 5.3 Functional felationships Between 
The State and Observation Vectors 


Obviously, some fairly involved chatn-rule Gititereneracion 


will ke required to ccmpute the partial derivatives in tne h 
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Matrix of equation 5.23 . Figure 5.3 summarizes the func- 
tional relationships ketween (kK, X, 8) SONG bee By the chain 


rule then 
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All that remains is to perform the reguired partial differ- 
entiations on equations 5.24 for substitution as appropriate 


Meme uatbrolS a.25) through 5.27 £0 that 
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It is left to the computer to calculate numeri- 
cally the partial derivatives described by eguations 5.28, 
and then substitute those values into equations 5.25 through 
5.27 to evaiuate the i matrix terms. Sequentially, these 
computations take place after the state and covariance esti- 
mates are calculated, and before the filter gains are 
computed for each iteration k, using % (k] k-1) in all casesz 


This ccmpletes the observer function derivations. 
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D. DISCUSSION OF EKF TRACKING KESULTS 
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From Range, Alpha, and Beta 


Bsing the DISSPLA graphics subroutines at the program's 
end, several sets of curves were generatec. Of general 
interest is the tracker's accuracy when R,X, and @ are all 
availiable, as opposed to R and K only, or X and Sonly 
(passive tracking, rcetentially). The series of graphs that 
comprise Figures 5.4 through 5.6 show the filter's perforn- 
ance for a single pass of the satellite. 

From these three graphs, 1t appears that the filter 
estimate does converge to the true target trajectory in all 
three cases, and that the filter eStimate is more accurate 
when all three variables are present than when only two are 


avaiiatkile. A single run of a random process haS a minimal 
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Figure 5.5 True and Estimated Target Track 
From Range and Alpha Only 


amount of statisticai Significarce, however; thus repeated 
trials of the tilter*s tracking will be Scequmrreda. in 
[Ref. 3], Bar-Shalom discusses in several instances the need 
for multiple sequential testing of any filter based on 
Stochastic principles; the statistical average of these 
"Monte Carlo runs" then begins to have some significance, 
increasing in value in direct correlation to tne number of 
runs considered. 

First, however, a criterion must be established for 
CrEOr CR vanrson. Bar-~Shalom's suggestion is the use of a 


State-sguared error function, defined for eack k as 
SSE(k) = KC) (KIRK) KO), 


Application of this criterion to €hiS “filter (note oe 


here) indicates that while the filter is able to track the 


rte Gere PA AY LY ee eS Nee aeSEC TCO 


NORTHERLY MILES 





= 
LEGEND 

° O- TRUE TPRGET TRACK 
5 Oo - ESTIMATED TARGET TRACK 
= 

S 

2 

S 

5 

2 

R | 


0 


t 

| 

| 

aa a am a a | 

1262.C 1253.0 1254.0 1265.9 i266.0 1267.0 1268.0 1269.0 1270.0 1271.0 1272 
spo} Pac ie eas 


Figure 5.6 True and Estimated Target Track 
Frog Alpha and Beta Only 


target in all cases, there is room for improvement in its 
perrormance. Specifically, the quality of the filter opera- 
tion can be gauged by how statistically well the SSE values 


ct 
Pits within a ‘confidence interval, where 
N = (number of runs) X (error vector dimension), 


In this case, when a 95% confidence interval (1.e€. the 
region within which 95% of the SSE values Should fail Lased 
on the pe Pamela ty GiStr@rucion) iS considered, less 
than half of the SSE values fall within it for any of the 
three cases under consideration. Speculation as to now tne 
filter's performance can be improved is reserved for Chapter 
Vi; ain the neantime, however, a Simplified criterion can be 
employed here to attempt a rough evaluation of the relative 


Cmrectiveness aft tne three Grier erent modes of filter 


operation. 
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Figuue,oa7 Simplified Error Criterion vs Time 
For a Single Trial of the Pilter 


Rh Simplified errer criterior that suggests itseif pight 

ke 
Sose(k) = ~~. (ae SG) nM 

wnere the time step functions as a scaling factor to roughly 
eguate tne significance of the fosition and velocity errors. 
Figure 5.7 shows the applicaticn Of Eh{S criterion eee 
Single rune 

Once again, the results seem to validate the hypothesis 
that the filter works better for all three variables than 
for two, but this is stihl ingorpmatitonmgericd ir onsen) ea 
Simo@dlet caval. Now that the error criterion has been estab- 
lished, however, tne Monte Carlo runs can be considered. 

The numbee of runs®is @hosen to be 20) and tie se 


application of the repeated trials was to smooth out the 


94 


mecarmwwywywu be rr 4 fet Wet FF jG FFF 


2710 


LEGEND 
O- TRUE TRRGET TRACK 
OF ES AS feo eee! “TRACK 


2706.0 2708.0 


NORTHERLY MILES 
2704.0 





2702.0 


! 

Te ee ee eee 

1282.9 i363.0 Loe .3 1285.0 1266.0 1267.0 1268.90 1269.0 1270.0 1273.0 1272. 
BeSncREY MiLES 


2700.0 


Figure 5.8 True an, Estimated Target Track 
FLOR venty" 10 ha, and Beta 
For Twen en € Carlo Runs 


estimated trajectories in Figures 5.4 through 5.6; the 
result was the series PoLogwbes 2.6 tnreugh 5.10 2 Ine 
eee rerceptions of the filter's ability to track the 
target appear to have been borne out, that is, the estimate 
does converge to the true track in the general case, and the 
filter tracxs better when all three variables of observation 
are available. However, the initial error in the estimated 
Petes )Phoved NOL § Lo be an anomary of the particular trial, 
but showed up repeatedly throughout the Monte Carlo runs. 
This suggestS an imbedded error in the initialization 


Eom Soeur O Ge thewt Liter OSs ite. me thewanrtial covariance 
’ p 


estimate; @uemewould aisSo help to explain the filter's 
fazpiwee tO place the SSE values ih the V,. confidence 
interval. Peamematron OL Figure 5.11, which shows’ the 
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Simplified error Criterion for vthe 20 eens, aiso tends to 
validate the results of the Single trial, and further 


Supports the supposition that the filter correctly tracks 
the target in each case, although there is some initial 
error. Finally, Figures 5.12 and 5.13 were drawn to take 
another closer look at the guality of the filter estiaates 
cver zQ runs. It now becomes cbvious that a primary source 
of €rror iS in the initial velceity estimates. Also, the 
filter's slow (30 time steps or 7.5 minutes) convergence to 
the true velocities indicates that Some werror inthe mnie 
covariance estimate is indeed fresent. Ironically, because 
one initial velocity estimate is low and the other is figh, 
the filter's initial estimate of the target's speed is 


fairly accurate; it 1S the e€Stiltlatevon tanget Ccounse me acer 
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Geeta lly in significant error. trtisecould well be @ func= 
tion of the geonetry of tne nerlinearities of the protien, 
indicating a potentiai pitfall cf assuming that white noise 
will remain white after undergjcing several nonlinear trans- 
forraticns. Mesemrng targets moving in other directions 
Brget wel, Shed lignt on the veracity of this supposition. 
Pomme 2emmmrerly, ditheugh the f11lter*s initializa- 
tion might be improved, it nonetheless demonstrates an 
ability to track the target in all three nodes of operation. 
It provides evidence that the taryet track is more accurate 
when all three variakles of okservation are available as 
opposed to only two, but that operating in the X- 4 mode 
does not really represent a Sigrificant degradation over the 


R-~ mode. This 1S an important result because it allows for 
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the pessibility of passive target tracking. The supposedly 
unobservable R-4 mode of operation is not attempted because 
the solution of the equations that comprise the filter algo- 
EPtnn Wiviout nh OF ye 1s difficult but feasible, but without 
using & it becomes impossible. Verification of the unob- 
Servability of the target tracked by range and altitude 
angle only thus requires a diiferent filter. 

This completes the discussicn of the results of the EKF 
tracking a surface target, and concludes the analytical 
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During the evolution of this thesis, it has become very 
clear why the observability of nonlinear systems has histcr- 
ically teen an unsclved problem, at least in the most 
general case, and inroads have reen made only gradually by a 
relatively few individuals. The algebra involved in 
applying any mathematical form of analysis to a nonlinear 
system is freguently so prohibitive as to make sich proce- 
dures of gueStionable value. Conseguently, the observ- 
ability analysis of a nonlinear system is frequentiy Lased 
as much on intuition as demonstrable mathematical j;rocf. 
Thus while the statement of this problem--observabiiity 
analysis for a satellite tracking an earth-bound target--is 
reasonably simple, it quickly becomes obvious its formula- 
tion Gs Mier. S peGisiaieal va every rererence available 
pertaining to either the sateilite-taryet dynamics or nonli- 
near ohservability uses a different notation, and there are 
very clearly defined crossover points in the information 
availiable. fost scurces either discuss the satellite's 
orbit in the language of physics, describe target observa- 
tions in a purely empirical sense, or discuss noniinear 
observability as applied to a general State srace. 
Conseguently, one of the major objectives of this thesis 
becomes the synthesis of aS much available informaticn as 
possible on the entire problea, and the Hormnulation, ota 
pertinent and coherent state vector model. For the sake of 
their universal recognition, Latitude and longitude are 
retained as the earth's surface coordinates, although much 
debate nas taken place on the subject of alternative coordi- 
nate svstems that might pertain. Certainly, one would hope 


that there exists another reference frane wherein the 
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mathematical observability anaiysis of Chapter IV could be 
Simplified such that the transfcrmation to rectangular coor- 
Ginates wouid not be required. SDoLLagG vee tne calculations 
required to itpiement the filter alyorithw of Chapter V 
might aiso be streamlined. Because there is potential for 
continued exploration of these and other aspects of this 
problem, the suggestions for further work made below are no 
less important than the specific conclusions that can be 


Obameetronm this research erfort. 


A. CONCLUSIONS 


AS 1S so frequently the case, the original intent of 
this research was overly ambitious in attempting to solve 
both tne complex nonlinear observability and target tracking 
probiems. The results presented here, nowever, do systemat- 
icaliy snow the interrelationship of taese problems, such as 
through the computation of the Jacobian matrices for the 
perturbation analysis, fOr the EKF,; ang 2Or the observ- 
ability analysis. This points cut the need for an efficient 
algorithm for high-order Jacobian evaluations. 

The analysis presented here forms a good base _ for 
continued research on tracking satellite deployment geometry 
and desired receiver information. The parametric descrip- 


tion cf the satellite dynamics pertinent to this problen 


involves a great deal of information not necessarily 
familiar to most engineers. This research proved inter- 
esting, and is well documented in the literature. There is 


awealth of diverse information available in sources such as 
[Ree- Sipand [Ref. 6 }. Tne intent here, however, 25 £6 
integrate this information intc a coherent overview of the 
mechanics of an orbitiny satellite and its observation of a 
taryet moviny on or near the Earth's surface, providing 
enough information to fully define the uroblem without 


drowning it in minutia. 
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Next the physical system is redefined in terms of the 
State vector space more familiar to Electrical Engineers. 
Also, the analysis establishes the computational procedure 
for transiation frcem observed satellite data to eéarth 
Ssucface coordinates. Appropriate coordinate selection here 
is obviously crucial to the ofservability analysis. ANY 
method of locating position cr describing motion on a 
Sphere, however, requires a radial vector and two angles, 
which is precisely what is defined by latitude, longitude, 
and altitude (i.e. the radial vector from the Earth's 
center). The further conversion from latitude and longitude 
to northeriy and easterly distance is made in an effort to 
at least partially linearize the problem; the merit of this 
decision 1S proven by Simulation, wnerein the constant 
velocity target is modeled very Simply in these coordinates. 
Had latiude and longitude been retained in the plant coordi- 
nates, even this simple case would have reyguired a highly 
nonlinear description. 

With regards to target observability, the Simpler cases 
considered first are intended to provide some intuition 
about the probiea. From the linearized case, it is deter- 
mined that having the position variables available for two 
or more cbServations is equivalent to observing both posi- 
tion and velocity Guring 2 simolertoon From the case cf 
the nonlinear model of the stationary target, it becomes 
obvious that the coordinates derived in tne earlier charters 
are not especiaily well-suited to this method of observ- 
ability analysis. Thus another transformation 1s performed, 
this time fron the spherical space into a rectanguiar one. 
In the process, it is discovered that a stationary target's 
position can be fixed with twe or more looks at a single 
variable of observaticn, which migitt Prove of jreat ulm ae 
if itis desired to tocGate passively “a Vand-tascdetacdan 


transmitter or communications station. When the moving, 


airborne, general-case target is ultimately considered, it 
1s shown that such targets prcve unobservable when only a 
Singie parameter of observaticn is available. When the 
variables cf opservation are considered in combination, 
however, the target becomes otservable in virtually every 
case that includes at least one position variable. The 
insight developed during the mathematical derivation of the 
coordinates suggests that the R-4 combinations should prove 
the exception to this rule, and lead to an unokservable 


system whether considering complete state or just position 


observability. Demonstrating unobservability conclusively 
by this method, however, would regquire the evaluation of 
some 130,000 determinants; thus applying tne analytical 


method to demonstrate this result is of questionable practi- 
Geli t y. While a complete nonlinear obsServability analysis 
then is ailimost impossible, the analytical method does 
provide further insight into relevant aspects of observ- 
ability for the probilen. One farticularly important result 
of the analysis is the c-only tracking of the maneuvering 
surface target, and the «4 tracking of the generai 
target, both very useful if a fassive sensor is of tactical 
advantage. As a footnote, the greater sensitivity or the 
oLserver matrices to fluctuaticns bg VS) than to Ror x iS 
Giscussed orierly; further development of this topic is 
certainiy possible, although time did not permit it here. 
The first-cut attempt at a tracking filter design serves 
to validate many of the concepts and intuitions brought to 
light during the observability analysis, micluding =» the 
feasibility of xX-f tracking, and the greater accuracy 
possible when employing all the variables o£ observation. 
The inherent problems of the EKF in general and this version 
Mme aeeweulate lead cne to conciuce that this 1s not the 
Optimum metnod for tracking a target. Fortunately, other 


methods are available, including the Truncated Seccnd Order 
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Filter (ISOF) and the Gaussian Second Order Filter (GSOF), 
both discussed in [ Ref. 9}. Their greater utility snould 
not, however, take away from the fact that this EKF does 
generate an estimate that converges to the target's actual 
track in all three phases of its operation, at least for the 
Smail number of cases considered here. It aiso proves that 
the coordinates developed in Chapters II and IIi to describe 
the proklem, while not particularly well-suited to the 
observability analysis of Chapter IV, work demonstrably well 


in engineering application. 


Be. RECOMMENDATIONS FOR FURTHER WORK 


During the development of this thesis, it became obvious 
that to discuss completely ail aspects of this problem at 
this time would not ce possible; thus these recommendations 
for further work are made. A final word on the subject of 
coordinate selection: for every eguation that appears in 
Appendix A as part of the Jacokian observer matrix calcula- 
tions, at least half again as many equations were derived 
during abortive attempts to apply this method of nenlinear 
observability to the problem in other coordinate systems, 
including latitude-longitude-alriride yy a) 2) ema R= X50 as 
the plant (i.e. system) coordinates. None of these provide 
any Simpler mathematical results; in fact, they all frove 
even more complicated to employ. Perhaps by redefining the 
problem in terms of a new set cf coordinates not discussed 
in this thesis the observability analysis could become nore 
Streamlined and/or definitively conclusive, but no such 
coordinate system has revealed itself thus far. 

Next the questicn of the tracking filter should Ge 
addressed. Had more time been availiable, this report would 
have attempted to improve the filter performance, either by 


eradicating the inherent errors apparent in the state vector 
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and covariance matrix initializations, or elise using a more 
sophisticated filter such as the TSOF or GSOF. Further, a 
two-stage noniinear filter has recently been proposed which 
ShowS superior tracking performance to any of these; its 
Specifications and performance are discussed in [Ref. 10] 
and [Ref. 12}. The employment of these or other more 
complicated filters might well afford the opportunity to 
verify the feasibiiity of an g@-only track of the maneu- 
vering surface target or the furported unobservability of 
bee R- Gebra che The goal in this thesis, however, is not 
@etmaliy to design the Optinut tracking filter, but to 
demonstrate the possibility of doing so. The EXF is a rela- 
tively simple algorithm, and thus is useful here; even so, 
there is room for further analysis of its performance. 
Other methcds of initialization might be considered, partic- 
ulariy for the covariance martix; PSO mC e fe tbacKkil ng OL 
other targets, such aS a maneuvering surface target cr an 
air target with its six-dimensicnal state vector, mignt be 
attempted. An analysis of its sensitivity could be 
conducted, varying bcth the system and the observer noises 
in turn; this would hopefully validate the declaration in 
Chapter IV that the system shouid be more sensitive to fluc- 
Cia nens ineyf tham those in R cr x. All of these varia- 
tions on the theme of filter oferation could of course asso 
be undertaken for the alternative filters. 

No mention has been made thus far of multiple targets or 
multiple observers, cut of course the situation explored in 
this thesis is the Simplest fessible, that is, a Sinyle 
target being observed by a single satellite with a perfect 
Prabapi vat (OL observation when in view. Pome tuality, oF 
course, the probiem is far more complicated, as the consid- 
eration of multiple targets and multiple observers immedi- 
ately introduces the probability of false target detection; 


of igcroper Gorrelation of targets to Successive or 
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Simultanecus observations; and of missing detection aito- 
gether. By viewing the entire process start-to-finish as 
stochastic rather than deterministic, the problem ultimately 
achieves a true real-world significance. In @§ Ref. a, 
Bar-Skhaiom has pulled together a wide variety of resources 
to anaiyze just this sort of situation in the yeneral sense; 
the application of the principles discussed there for this 
particular problem would arguably generate several more 
reports aS involved as this one. 

On a higher level, the theoretical implications of the 
applicability (or limitations thereof) of the methcd of 
fRef. 1) to this preblem and others of equal difficuity 
Hight be discussed. Previous é€mployment of this method has 
Shown the utility of such observability testing in deter- 
Mining the lack of observability of the bearing-only target 
(BOT) for typical earth-borne radar and sonar applications. 
The compiicated geometry of this problem has not yet yielded 
such complete conclusions; the potential exists, however. 
Another topic that might be exfiored further is the verva- 
Sive nature of the Jacobian matrix, which appears as part of 
the observability analysis; in the development of the nonli- 
near observer here and elsewhere, i.e. [Ref. 10], 
[Rei. 11], and [Be£f. 12], and again in the computaticn of 
the observer functions for the filter. Finally, tne contin- 
uing effort must be made to resolve all these results, 
current and future, with the empirical evidence available of 
track observability and maintainability, taking into account 
both tne capabilities of the sensor and the wiles of the 


target. 
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DERIVATION OF OBSERVATION MATRICES 


The foliowiny series of derivations follows the pattern 
established in Chapter 4 Section B for establishing the 


terms in the Jacobian observer ratrices for the general case 


tanget. The accelerationS are assumed to be piecewise 
Constant (1. ¢. constant for each time step of a discrete 
observer); thus third and higher order derivatives of the 
position variables will be zero. The derivations are 


Geelped by two's into Rk and ay > and ue , and 4 and 4 ONE 
the reason that nany of the terms in their natrices are the 


salle. 


A. OBSERVER MATRICES FOR RANGE AND RANGE RATE 
@omeeder first that. in the rectangular coordinate 
SYStem, 
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Taking the tartial derivative of each of the above equations 


with res ,;ect to n, 
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Taking the partial derivatives c£ tne first set of equations 


Bith resrect to n, 
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Somme LOL the nighest order partial derivatives as before. 
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ard n; combining those results will lead to the first 


observer matrix de Ge the R observation: 
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last equation of each set of Seven partial 


derivatives was not required at this time; they will become 


necessary in tne case of R. The lower rows of this matrix 
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The resulting matrix is 


nN oe oho 
navn ee hh 

2 Cte mete Gu Gh 
en Omen ace OC CH 
Ge 2 Oo 0 


\! 


Oe NGmnO=s Oxe©- © 


- 


mee ene “eres sor 6the R JWatrix J, mave already been 
determined; they come from the same sets of equations just 
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Be. OBSERVER MATRICES FOR AZIMUTHAL BEARING AND 


BEARING RATE 


The derivation of the Jacobian observer 


AZIMUTHAL 


matrices for the 


azimuthal bearing variables ™ and X& follows the method of 


the previous section, except that the result is not guite so 


algebraically nice. Consider first the definition of ips} 


this coordinate system: 
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Gakingetne partial deriyative of each of the first seven 


e,uations with respect to n follows the same pattern; 


solviny again for the highest order derivative, 
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Obviously, all the partial derivatives with respect toh 


will ke zero. 
This leads to the formation of the J, and J, matrices 


X, A de dA 9 0 yi wht el, He ONE 


Je [Are Ad OO}, | KK Lk co 
a Cn oe ae a 


CS 
\N Ww \Nv iN 
Cen ha CO vo ovo 8 
ok i Kr OA, L266 
on oe ol € O Cy “A ~ iy ONS 


which are simplified Fy the same sort of row-reducing method 


Sapeoyedein the J, and J_ matrices. This time, however, all 
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C. OBSERVER MATRICES FOR ALTITGDE BEARING AND ALTITUDE 
BEARING RATE 


The derivation or the J. and Jo matrices for tha alti- 
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PORTE SAPGELIZE-TARGET SIMULATION AWD THE 
EKF ITMNPLEAENTATION 
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